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SUMMARY 


The  analytic  modeling  of  well-organized  rotating  convective  storms  is 
invoked  to  elucidate  the  evolution  from  a  moderately  intense  one-cell 
vortex,  characterized  by  low-level  pressure  deficits  on  the  order  of  one 
percent  of  atmospheric  pressure,  to  a  ve"y  intense  two-cell  vortex,  char¬ 
acterized  by  low-level  pressure  deficits  on  the  order  of  ten  percent  of 
atmospheric  pressure.  The  physical  distinction  between  the  two  stages 
rests  on  the  insertion  of  a  dry,  compressionally  heated,  nonrotating, 
central  downdraft  of  originally  tropopause-level  air  in  the  transition 
to  the  more  severe  case.  The  quasisteady  mature  description  of  the 
thermohydrodynamic  structure  of  each  vortex  is  being  developed,  and  then 
the  conditions  for  transition  from  the  moderately  intense  to  the  very 
intense  vortex  will  be  sought.  The  practical  motivation  is  to  make 
progress  toward  the  highly  desirable,  but  very  formidable  task  of  being 
able  to  anticipate  which  tropical  storms  or  minimal  hurricanes  will 
evolve  to  supertyphoons . 

With  a  description  of  the  one-cell  storm,  believed  adequate  for 
present  purposes,  already  available,  description  of  the  two-cell  vortex 
is  the  current  challenge.  In  particular,  since  the  potential-vortex  and 
surface-inflow  subdivisions  of  the  structure  of  a  very  intense  vortex 
are  in  hand,  a  realistic  picture  of  the  updraft  and  its  location  are 
sought . 

The  modeling  proceeds  from  basic  thermohydrodynamic  principles. 
Nonessential  geometric  detail,  as  well  as  association  of  conclusions  with 
the  details  of  particular  parameter! zations ,  is  being  avoided  as  much  as 
possible. 

Tentative  results  achieved  thus  far,  by  an  integral-type  treatment 
that  seeks  the  location  of  the  eye  wall  and  its  properties  with  height, 
point  to  a  particularly  sensitive  dependence  on  the  vertical  variation 
of  the  angular  momentum  in  the  portion  of  the  vortex  outside  the  eye 
wall.  Perhaps  for  only  certain,  narrowly  defined  stratifications  of 
angular  momentum  may  the  conservation  laws  be  compatible  with  two-cell 
structure  for  an  intense  atmospheric  vortex. 
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INTRODUCTION 


1  . 

I . 1  Statement  of  the  Objective 

The  principal  contribution  still  to  be  made  in  forecasting  the 
intensity  of  tropical  storms  is  to  identify  the  observable(s )  that  permit 
one  to  anticipate  (as  much  in  advance  as  possible)  which  incipient  storm 

M 

will  become  hurricanes  and  which  will  not. 

Very  roughly  about  half  of  all  tropical  stroms,  one-cell  vortices 
extending  vertically  from  sea  level  to  tropopause  and  radially  many 
hundreds  of  miles,  transform  into  two-cell  vortices  (hurricanes)  (Fendell 
1974).  The  transition  involves  the  insertion  of  a  column  of  relatively 
dry,  clear,  nonrotating,  tropopause-1 evel  air  at  the  center  of  the  vortex. 
Accordingly,  the  region  of  intensely  swirling,  very  cloudy,  torrentially 
raining  updraft  becomes  an  annulus  surrounding  that  "eye"  of  dry  air; 
i.e.,  an  "eye"  has  then  been  inserted  within  an  "eye  wall".  A  formidable 
vortex,  with  a  peak  swirl  of  less  than  100  mph  and  horizontal  pressure  . 
deficit  at  sea  level  of  roughly  twenty  millibars  (the  tropical  storm,  as 
defined  here),  can  be  transformed  into  a  terrifying  vortex  with  a  peak 
swirl  of  some  200  mph  and  a  pressure  deficit  of  appreciably  more  than 
50  mb  (the  hurricane,  as  defined  here).  An  important  question  is;  "what 
telltale  observable,  accessible  to  aircraft,  satellites,  etc.,  would  per¬ 
mit  one  to  anticipate  the  transition  from  tropical  storm  to  hurricane?" 

The  intensity  of  a  tropical  cyclone  is  not,  of  course,  a  simple  monotonic 
increase  with  time  to  a  peak,  following  by  a  simple  monotonic  decay; 
rather,  there  may  be  several  transitions  back  and  forth  between  one-  and 
two-cell  structure,  suggestive  of  partial  insertion,  then  removal,  of  an 
"eye"  (Carrier  1971a). 


#Gray  (1979,  p.  51)  comments:  "Tropical  cyclone  forecasters  generally  agree 
that  skill  at  operational  forecasting  of  tropical  cyclone  intensity  change 
is  nearly  zero."  He  notes  that  trying  to  relate  intensity  to  cloudiness 
configuration  (as  observed  by  passive  satellites)  is  subject  to  signifi¬ 
cantly  erroneous  results,  and  furnishes  inadequate  detail.  It  may  be 
noted  that  the  gustiness  of  winds  is  not  discussed,  since  treatment  of  the 
variance  seems  premature  when  prediction  of  the  mean  is  still  not  possible. 
Gustiness  factors  for  tropical  cyclones  (ratio  of  peak  transient  wind  to 
wind  sustained  over  a  scale  of  one  minute  to  five  minutes,  depending  on 
the  correlator)  range  from  1.0  to  1.8  in  the  literature  (Brand  et  al  1979). 
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Normally,  quantifying  any  details  about  structure  in  tropical  cyclones 
is  highly  uncertain  because  proper  formulation  of  turbulent  diffusion, 
cummulus  convection,  and  radiative  transfer  is  unknown.  These  three 
important  phenomena  may  well  not  have  satisfactory  representation  in 
terms  of  macroscale  variables,  so  the  attempt  at  such  parameterization 
incurs  uncertainty.  Yet  such  parameterization  is  entailed  in  virtually 
all  contemporary  tropical  cyclone  modeling. 

1 .2  An  Approach  to  Two-Cell  Structure  for  an  Intense  Atmosphere  Vortex 

Whereas  most  contemporary  hurricane  modeling  entails  direct  numerical 
assault  on  all-inclusive  mathematical  formulation,  here  an  alternative, 
approximate  analytic  treatment  is  utilized  as  much  as  possible  (Carrier 
1970;  Carrier,  Hammond  and  George  1971).  The  hurricane  is  subdivided 
logically  into  those  parts  in  which  different  physical  phenomena  dominate; 
when  a  particular  subdivision  is  treated,  only  that  subset  of  the  full 
equations  that  retains  the  locally  dominant  physical  processes  is  solved. 
Such  a  subset  is  usually  more  tractable  to  solve,  and  parametric  variation 
easier  to  examine.  A  composite  global  solution  is  synthesized  from  solu¬ 
tions  to  the  subdomains  by  demanding  appropriate  continuity  of  dependent 
variables  and  their  fluxes  at  interfaces. 

Under  the  high-speed  portion  (only)  of  a  severe  rotating  storm,  the 
nonlinear  "surface"  inflow  layer  is  inviscidly  controlled;  only  in  a  very 
small  sublayer  (immediately  continguous  to  the  ground)  of  the  inflow  layer 
does  diffusion  enforce  the  no-slip  boundary  condition  (Carrier  1971b; 
Burggraf,  Stewartson  and  Belcher  1971;  Carrier  and  Fendell  1978).  The 
diffusive  sublayer  becomes  thinner,  and  the  outer  inflow  layer  thicker, 
as  one  moves  toward  the  center  of  the  vortex.  [An  analogous  result  is 
familiar  in  nonrotating  contexts  to  aerodynami cists  in  the  form  of  the 
following  bit  of  well-known  empiricism  (e.g.,  Launder  1964);  a  near-wall 
diffusive  layer  becomes  thinner  in  the  direction  of  a  favorable  (i.e., 
accelerating)  pressure  gradient.]  At  sea  level,  the  pressure  is  high 
at  the  lateral  edge  of  the  intense  vortex  and  the  pressure  is  greatly 
reduced  at  the  center,  such  that  a  strong  positive  pressure  gradient 
exists;  thus,  the  low-level  swirling  flux  (that  erupts  up  the  eye  wall) 
is  predominantly  inviscidly  controlled  as  it  proceeds  toward  the  hurri¬ 
cane  center.  The  cited  pressure  gradient  is  radial,  so  the  statement 
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of  a  limited  role  for  frictional  effects  pertains  primarily  to  the  state¬ 
ment  of  conservation  of  radial  momentum.  However,  the  azimuthal-velocity 
profile  is  largely  developed,  such  that  further  evolution  of  the  surface- 
inflow-layer  profile  near  the  center  of  the  vortex  is  not  to  be  anticipated. 

Radiative  cooling  is  important  only  in  the  outer  and  less  cloudy 
portions  of  the  storm  (Fendell  1974),  and  the  simple  well-defined  model 
of  moist-adiabatic  ascent*  suffices  to  describe  the  locus  of  thermodynamic 
states  in  the  eye  wall.  In  the  previous  paragraph  it  is  pointed  out  that 
the  low-level  flow  entering  the  eye  wall  is  nondiffusive.  Two  of  the 
three  processes  that  entail  parameterization  (specifically,  turbulent 
diffusion  and  radiative  transfer)  do  not  enter  significantly  in  the  por¬ 
tion  of  the  hurricane  flow  field  critical  to  quantitative  analysis  of 
two-cell  structure,  and  the  third  process  (cumulus  convection)  can  be 
formulated  locally  in  a  well-accepted,  well-defined  form. 

The  upshot  is  the  following  good  fortune:  one  of  the  few  questions 
about  hurricane  structure  that  can  be  definitively  formulated  by  con¬ 
temporary  modeling  happens  to  be  one  of  the  most  important.  Specifically, 
that  question  is,  what  is  the  balance  of  forces  and  thermodynamic  states 
in  the  central  region  of  a  tropical  storm  that  permits  two-cell  structure 
to  arise  as  an  alternative  to  one-cell  structure? 

1 . 3  The  Boundary- Val ue  Problem 

For  axi symmetric  inviscid  swirling  flow  in  which  a  radial  influx 
becomes  an  axial  upflux,  the  boundary-value  problem  is  well-defined. 

If  one  confines  attention  at  the  outset  to  low  altitudes  such  that 


★ 

Water-vapor-contai ni ng  air  rises,  expands  and  cools  dry-adiabati cally 
(at  constant  water-vapor/ai r  mixture  ratio)  until  (in  the  absence  of 
supercooling)  the  temperature  drops  far  enough  to  imply  saturation. 
Thenceforth,  it  continues  to  ascend  such  that  saturation  of  water  vapor 
is  maintained  at  the  local  thermodynamic  state,  the  condensed  excess 
water  substance  falling  out  so  that  is  is  unavailable  for  re-evaporative 
cooling.  (Hence  the  process  cannot  be  reversed.)  The  originally  low- 
level,  convectively  unstable  fluid  undergoing  this  ascent  is  taken  to 
rise  too  quickly  to  entrain  any  ambient  air,  and  hence  remains  unmixed 
with  cooler  or  drier  air.  In  fact,  of  course,  some  mixing  with  nearby 
ambient  air  occurs  and  some  re-evaporation  into  this  unsaturated  mix  also 
occurs,  but  the  no-mixing  case  is  realistic  enough  to  provide  valuable 
insight. 
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incompressibility  suffices,  the  formulation  is  available  in  texts  (e.g., 
Batchelor  1967).  Angular  momentum  is  constant  on  streamlines,  continuity 
is  enforced,  advective  and  centripetal  accelerations  balance  the  radial 
pressure  gradient,  and  the  total  pressure  head  is  constant  across  stream¬ 
lines  {Carrier,  Dergarabedian  and  Fendell  1969).  The  annulus  is  bounded 
on  one  side  by  a  free  streamsurface  that  demarcates  the  eye  wall  interface 
with  the  eye,  and  on  the  other  side  by  a  free  streamsurface  that  demar¬ 
cates  the  eye  wall  interface  with  the  outer  vortex.  The  position  of  both 
these  bounding  streamsurfaces  must  be  found  in  the  course  of  solution, 
under  enforcement  of  statements  concerning  the  continuity  of  pressure. 

One  result  sought  is  the  displacement  of  the  eye  wall  from  the  axis  of 

symmetry,  which  is  also  the  axis  of  rotation.  It  is  anticipated  that 

% 

this  displacement  ultimately  increases  with  height.  The  mathematically 
elliptic  character  of  the  boundary-value  problem  reflects  the  physical 
fact  that  streamwise  and  transverse  derivatives,  and  streamwise  and  trans¬ 
verse  velocity  components,  are  of  comparable  magnitude;  there  is  no  single 
preferred  direction  of  change. 

Useful  insight  is  attainable  at  the  sacrifice  of  detail  by  using  an 
"integral  method"  (Finlayson  1972),  or  one-dimensional  approach.  Roughly, 
one  "forces"  a  parabolic-like  character  on  the  problem  by  effectively 
integrating  out  transverse  dependence  and  then  seeking  streamwise  variance 
from  the  simplified  boundary-value  problem. 


The  increasing  displacement  of  the  eye  wall  from  the  axis  with  altitude 
seems  likely  under  conservation  of  angular  momentum,  since  density 
decreases  with  altitude.  Conservation  of  mass  also  suggests  that  the 
cross-sectional  area  of  the  eye-wall  annulus  increases  with  altitude. 

The  fact  that  the  eye  wall  tilts  radially  outward  with  increasing  height 
is  suggested  further  by  certain  satellite  photographs  of  hurricanes 
(Fendell  1974).  If  the  eye  wall  were  vertical,  pressure  deficits  achiev¬ 
able  under  compressional  heating  in  the  eye  could  not  serve  as  an  explana¬ 
tion  of  intense  swirl  speeds  observed  in  the  eye  wall  (cf.  Malkus  and 
Riehl  1960). 
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2.  A  BOUNDARY-VALUE  PROBLEM  FOR  THE  TURNAROUND 

An  inviscid  incompressible  steady  axi symmetri c  model  of  the  turnaround 
region  in  a  severe  vortical  storm  is  adopted,  in  cylindrical  polar  coordi¬ 
nates  with  origin  at  the  ground  on  the  axis  of  rotation  (and  of  symmetry) 
(Figure  1 ) . 

The  azimuthal,  radial,  and  axial  velocity  components  (v*,  u*,  and 
w*,  respectively),  in  view  of  conservation  of  angular  momentum  and  of 
conservation  of  mass,  may  be  written  [r*  =  r*V*,  where  r*,  is  the  (given) 
radius  of  peak  swirl  in  region  I,  and  where  V*  is  the  (given)  value  of 
that  peak  swirl] 


(r*v*) 2  =  r*2F(’^)  ; 

(2.1) 

_r*u*  =  rw  M.  . 

r  u  IT  3z*’ 

(2.2) 

*  *  r*h* 
r*w*  =  - . 

tt  3r* 

(2.3) 

Super  asterisk  denotes  a  dimensional  quantity;  no  asterisk,  a  dimensionless 
quantity.  The  cylindrical  radial  coordinate  is  r*;  the  axial  coordinate, 
z*;  the  maximum  height  of  the  surface  inflow  layer,  h*;  the  streamfunction 
for  the  secondary  flow  (involving  velocity  components  u*,  w*,  only) ,  ip; 
the  function  giving  distribution  of  angular  momentum  with  streamfunction, 

F. 

Conservation  of  radial  momentum  is  (p^  is  density) 


u* 


3u* 

3r* 


+  w* 


(2.4) 


Bernoulli's  equation  is 


+  w*2  .  2l  .  r*2F(») 


2r*2 


g*z* 


Pi 


2r  f 


(2.5) 
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in  (2.5)  the  Bernoulli  constant  is  taken  as  universal  for  all  streamlines 
emerging  from  region  II  in  view  of  the  known  relation 


[u*(r*,z*)]2 

_ 


+ 


r*2FU) 


(2.6) 


and  the  fact  that  the  pressure  above  hydrostatic  is  axially  invariant  for 
z*  <_  h*,  within  the  conventional  boundary-layer  approximation.  Within 
the  cyclostrophic  approximation, 

p*  -  P*  =  (p*/2)V*2  =  (p*/2)(r*2/rf),  (2.7) 

where  p*  is  the  qround-level  ambient  pressure  (given). 

3 

From  a  crude  treatment  of  the  turnaround,  to  be  published  elsewhere, 
(h*/r*)  =  0.2  for  relevant  values  of  r*,  V*,  where  typical  (given)  values, 

^  rn  a 

for  a  severe  tornado,  are  r*  =  160  m,  V*  =  100  m/s.'  In  any  case,  h*  and 
r*  are  taken  as  known.  Thus,  in  addition  to  its  use  in  (2.2)  and  (2.3), 
h*  is  employed  to  nortdimens  ianal  ize  the  independent  variables: 


x  =  ~r*/h*,  y  =  Tiz*/h*  . 


(2.8) 


Subtraction  of  the  radial  derivative  of  (2.5)  from  (2.4)  gives 


x 


+  4', 


yy 


(2.9) 


where  subscripts  x  and  y  denote  partial  differentiation.  The  swirl  profile 
F(v)  is  known  from  conditions  holding  at  x  =  xx  =  irr*/h*  =  u/h.  It  is 
known  from  solution  of  region  II  (Carrier  and  Fendell  1978)  that  F  is  a 
monotonical ly  increasing  function  of  ty,  such  that  F(0)  =  0  and  F(I)  =  1, 
where  the  datum  ip( x  >  xl50)  =  0  is  adopted,  and  from  (2.2)  and  (2.8), 


#For  a  hurricane,  (h*/r*)  =  0.1,  r*  =  15  km,  V*  =  100  m/s. 

i  i 
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(2.10) 


TT 

»tt)  =  /-(x, , 


C)d;  =  I  , 


where  u  =  u*/V*,  and,  for  completeness,  v  =  v*/V*.  Clearly,  for  F  linear 
in  ip,  (2.9)  becomes  Poisson's  equation,  and  for  F  quadratic  in  \p,  (2.9) 
becomes  Helmholtz’s  equation.  Although,  from  known  results  for  region  II, 
such  simple  forms  for  F(tf>)  are  not  detailed  replications  for  (r*v*)/F*  of 
interest,  still  their  tractability  urges  their  adoption. 

The  pressure  along  the  streamline  that  separates  from  y  =  0  at  x  =  xj 
remains  at  p*,  for  the  model  of  a  nonrotating  "eye"  isobaric  over  an  altitude 
range  in  which  gravity  implies  negligible  pressure  changes.  Then,  from 
(2.5),  since  F(0)  =0,  . 


<»„>’  *  <V2  s  (i) 


on  tp  =  0. 


(2.11) 


The  pressure  along  the  streamsurface  ijj  =  I ,  which  passes  through  the 
circle  x  *  Xj ,  y  =  -,  should  be  consistent  with  the  pressure  in  the  potential- 
vortex  region;  hence 


p(x,y)  *  p;  + 


(2.12) 


from  the  cyclostrophic  balance,  where  the  nondimensional ization 


(2.13) 


has  been  adopted.  From  (2.5),  (2.12),  and  (2.13),  with  neglect  of  the 
gravitational  term,  since  F ( I )  3  1, 


ij>2  +  y2  =  0  on  =  I . 
x  y 


(2.14) 


This  condition  merely  states  that  there  is  no  contribution  from  the 
secondary-flow  velocity  components  u*,  w*  to  the  pressure  field,  where 
region  III  interfaces  with  region  I. 

Equations  (2.9),  (2.11),  and  (2.14)  constitute  the  boundary-value 
problem  of  interest,  when  supplemented  by  a  relation  for  F ( ^ ) ,  obtained 
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from  previously  executed  analysis  of  region  II.  By  (2.1),  (2.2),  (2.6), 
and  the  fact  that  wfx^y)  =  0,  the  relation  for  F(^)  gives  implicitly  an 
expression  for  0  <.  y  <.  it,  with  which  to  initiate  the  analysis 

(see  below).  Finally,  it  is  anticipated  that  \jj(x,y),  as  given  by  the 
boundary-value  problem  is  periodic  in  y;  this  periodicity  reflects  the 
fact  that  the  formulation  must  be  revised  after  completion  of  about  one 
cycle.  More  explicitly,  it  is  anticipated  that  the  locus  of  any  stream¬ 
line,  say  ip  =  1/2,  as  y  increases,  decreases  to  a  minimum  value  of  x,  then 
increases  in  x  to  recover  its  initial  value  x  =  k1  and  to  achieve  a  peak 
value  in  x,  before  decreasing  back  to  x  =  X;  to  begin  another  period.  It 
is  to  be  expected  that  the  undulating  flow  so  generated  will  be  unstable 
somewhere  along  its  trajectory,  but  the  immediate  goal  is  to  find  that 
flow  and  discussion  of  the  stability  problem  is  deferred  until  that  is 
done.  Of  particular  interest  are  the  amount  by  which  v  exceeds  the  value 
associated  with  x  =  xj ,  o  -  I.  and  the  values  of  x  and  tp  at  which  the 
maximum  occurs . 

For  explicitness,  a  rough  characterization  of  the  function  ^(Xi.y), 

0  <  y  <  ^,  is  now  presented.  In  the  largely  inviscid  portion  of  region  II, 
(2.6)  may  be  rewritten  as 

u2(x,y)  +  v2(x,y)  =  (x^x)2.  (2.15) 

From  (2.1),  (2.2),  and  (2.15),  one  obtains 

O'Hxi  ,y)  (  r  1 2 1  r 

~u(xx ,y)  =  - -  =  j  1  -  [v(xi  ,y)J  j  =  j^l  -  F(^)J  2,  (2.16) 

where  0  £  F(^)  <1  for  0  <  ip  <_  I.  Thus,  i^(xj,y)  is  given  by 

*i  .y) 

dy 

jl  -  [F(m)]| 2 

where  ijj  -*  I  as  y  "  [cf.  (2.10)].  For  F(^)  =  W I, 

'P  =  y  (]  ‘  $n)  ’  1  =  f  =^v2(xl  =  ^  (1  ‘  2^)  •  (2-18) 


(2.17) 
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For  F(<j/)  =  (ip/I  )2  , 


V  =  2^sin(y/2),  I  =  2"2  =^v2(x1  ,y)  =  sin2(y/2).  (2.19) 


In  fact, 


v(xi,y)  =1  -  exp  (-  £  y )  ,  c  =  3, 

(2.20) 

might  be  more  realistic,  but  is  far  less  tractable. 

In  summary,  if 

V(x,y)  ,  G'W 

(2.21) 

then  the  boundary-value  problem  is,  for  0  <_  'F  <_  1 , 

*  (  7  \  )  *  \y  + 

(2.22) 

H'2  +  Y2  =  { -T— —  I  on  Y  =  0,  y2  +  Y2  =  0  on  Y  =  1 ; 

x  y  \ixi/  x  y 

(2.23a) 

'f(xi.y)  given,  with  periodic  in  y. 

(2.23b) 

One  might  consider  interchanging  the  roles  of  the  dependent  variable 
i|j  and  the  independent  variable  x,  and  using  perhaps  some  technique  within 
the  framework  of  the  method  of  weighted  residuals  to  extract  the  desired 
information  from  (2.21 )-(2.23) ,  However,  here  alternative  procedures  are 
used  to  study  (2.21 )-(2.23) . 
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3.  TREATMENT  OF  A  LOW-LEVEL  SEPARATED  LAYER  WITHOUT  STRUCTURE 


The  momentum  balance  in  an  axisymmetric  separated  boundary  layer 
without  structure  is  examined.  The  thin  sheed  (or  "eye  wall")  is  the 
demarcation  between  an  "eye",  isobaric  at  pressure  p*,  and  a  potential 
vortex  with  radially  dependent  pressure  given  by  (2.12).  At  r*  =  r*, 
p*  in  the  sheet  equals  that  in  the  "eye"  (Figure  2). 

Henceforth  in  this  section  the  symbol  r*(z*)  denotes  the  outside 
surface  of  the  "eye  wall"  without  structure.  At  the  point  A,  given  by 
[r*(z*),z*],  the  principal  radii  of  curvature  are  {1  +  [r*‘ (z*)]2}3/2/r*"(z*) 
and  r*(z*){l  +  [r*‘ (z*)]2 p/2 .  Super  prime  denotes  ordinary  derivative 
with  respect  to  the  argument  of  the  function.  The  velocity  component  in 
a  plane  containing  the  axis  of  symmetry  is  denoted  q*,  while  the  velocity 
component  (swirl)  in  a  plane  perpendicular  to  the  axis  of  symmetry  is 
denoted  v*.  Thus,  if  t  is  a  unit  vector  in  the  plane  containing  the  axis, 
and  $  is  perpendicular  to  ^  and  refers  to  the  azimuthal  component  in  a 
cylindrical-polar-coordinate  system,  then 

v*  =  q*^  *  v*§.  (3.1 ) 


The  force  balance  perpendicular  to  the  thin  sheet  at  [r*(z*),z*] 
equates  the  pressure  gradient  consistent  with  a  potential  vortex  outside 
the  sheet,  to  the  component  of  acceleration  perpendicular  to  the  sheet. 
There  are  two  contributions  to  the  acceleration,  each  involving  the 
square  of  a  velocity  component  over  an  appropriate  radius  of  curvature: 
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3h* 


q*2r*" 


,*2 


[l  +  (r*')2  3/2  r*[l  +  (r*1)2]1 


/  2 


(3.2) 


where  h*  represents  a  coordinate  running  across  the  sheet. 

Variations  in  velocity  occurring  across  the  sheet  are  not  resolved 
(T*  =  r*V*): 
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'  v*2dh*  =  r*2  /  - — - =  r*2V*2  /  - — - 

J  [r*(h*)]2  1  J  [r*(h*)]2 


r*2V*2A* 


r*2  ^  [r*(h*)]2 


(3.3) 


3*2dh*  =  V*2 


f[^]‘ 


dh*  =  V*2B* 


dh*  . 


(3.4) 


The  potential -vortex  form  of  v*  is  used  in  the  definition  of  A*. 

From  these  definitions  and  from  (2.12),  to  an  accuracy  that  serves 
current  purposes. 


0-S)- 


A*y*2r*2 


B*y*2r*n 


>[l  +  (r** )2] 1/2  [l  +  (r*1 )2J 


.  (3.5) 


The  following  nondimensionlization  is  introduced: 


x=H  7  =  11  a  =  lA*  b  =  2B^ 
i  l  ri  l 


(3.6) 


If  A*  =  (h*/2),  B*  =  (h*/2),  seemingly  reasonable  values,  then 


•  h*  •  h* 

ct  =  pr  i  h,  8  =  5*  =  h 


(3.7) 


Under  (3.6),  (3.5)  becomes  (a, 8  specified) 
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1 


1 


Ot 


8x" 


Z2  _  „  1/2  _  „  3/2 

X  X3(1  +  X'2}  (1  +  X'2) 


;”=i  (3-8) 


In  this  translationally  invariant  equation,  it  is  taken  that  x(z  =  0)  =  1. 
It  is  shown  below  that  x  is  periodic  in  z.  Solution  is  sought  for  positive 
and  negative  z,  where  the  boundary  conditions  are 

x(0)  =  1  (3.9) 

x 1 ( 0)  =  x-,  given  const.  >  0.  (3.10) 

The  boundary  conditions  preclude  odd  or  even  solution  for  x  and  z.  Sought 
are  x+,  the  largest  value  of  x,  which  occurs  where  x'  =  0,‘and  x  ,  the 
smallest  value  of  x,  which  also  occurs  where  x1  =  0.  In  that  x  <1, 
there  is  "overshoot"  of  the  swirl,  and,  from  (2.12),  there  is  decrease  of 
pressure  from  the  value  that  holds  in  the  "eye",  i.e.,  in  0  <  r*  <  r*(z*h 
it  is  reiterated  that  the  magnitude  of  the  :wirl  overshoot  and  pressure 
decrease  are  of  particular  interest. 

While  numerical  integration  of  (3.8)-(3.10)  is  required  ultimately, 
some  preliminary  treatment  is  helpful.  Since  (3.8)  is  translationally 
invariant  in  z,  phase-plane  analysis  is  introduced  (Figure  3): 


Thus , 


dx 

p<x)*i£i- 

dP(x) 

dx  dP 

P^ 

dz 

dz2 

dz 

dz  dx 

dx 

BPP' 

a 

+  1 

-  1  . 

-  3/2  ,  -  1/2 

(1  +  P2)  x3(l  +  P2) 


(3.11) 


(3.12) 
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The  slope  is  infinite  at  P  =  0  and  is  zero  where 


,  -  1/2 
X3(1  +  P2) 


+  --1=0 

x2 


(3.13) 


the  intersection  of  the  curve  at  (3.13)  with  the  P  =  0  ray  is  given  by 


a 


1=0, 


(3.14) 


where  x*  =  1  for  a  =  0,  x*  >  1  for  a  >  0.  For  1  >>  a  >  0, 

X*  =  1  +  (a/2)  -  (3a2/8)  +  ...  . 

Although  the  following  development  is  not  pursued  to  the  extent  of 
obtaining  results,  it  may  be  worth  noting  that  if 

H(x)  =  (1  +  P2)"1/2,  (3.15) 


then  (3.12)  becomes 


-SH‘ 


(3.16) 


If 


(3.17) 


then 

26  4!^  -  aH  =  x-l/2  -  I’3/2, 


(3.18) 


(3.19) 


where  E  is  a  const,  of  integration.  From  (3.15),  (3.17),  and  the  defini¬ 
tion  P  =  (dx/dz),  one  may  write  formally 
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dx 

dz 


1  +  [H(l)-2] 


1/2=  d( t~ 1 /2  ) 
dz 


(3.20) 


Thus , 


dz  = 


d( t*1/2 ) 


1  +  CH(t)]"2 


(3.21) 


where  H(t)  is  given  by  (3.19). 

For  the  special  case  a  =  0,  i.e.,  no  swirl,  multiplication  of  (3.12) 
by  x'  yields 


6 


0 


+ 


(x  -  D2  .  6 

-  m  ’ 


(3.22) 


where  the  constant  of  integration  has  been  written  as  (B/m),  with 

-  „ 

m  =  (1  +  x^2) 


(3.23) 


for  consistency.  It  is  recalled  that  x^(h  dx(0)/dz)  is  a  given  positive 

finite  constant.  Inspection  of  (3.22)  reveals  that  x'  is  maximum  at 
•1/2  -1/2 
x  =  1,  so  (1  +  Xq2)  is  the  maximum  value  of  (1  +  x'2)  '  ,  whence  the 

symbol  m.  At  x'  =  0, 


-+ 

X“ 


(3.24) 


For  x^  0  so  m  -*■  1 ,  x-  merge  to  unity;  this  case  involves  a  vertically 
separating  surface  inflow  layer,  and  hence  no  overshoot.  For  x^  »  so 
m  -*■  oo,  x-  remain  bounded; 


-+ 

X" 


(3.25) 


this  case  involves  effectively  horizontal  inflow  of  the  separating  surface 
inflow  layer,  and  leads  to  the  minimum  value  of  x  for  fixed  8.  For 
8  =  0.1,  x+  =  1.37  and  x  =  0.73;  for  8  =  0.2,  x+  =  1.56  and  x  =  0.64. 
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The  plausible  range  of  S  =  0(0.2),  from  a  crude  analysis  of  the  turnaround 
to  be  published  separately;  it  is  recalled  that  0  is  associated  with  in¬ 
plane  motion.  These  results  suggest  what  proves  to  be  a  general  trend: 
x  decreases  as  and  6  increase. 

Finite  values  for  a  indicate  finite  swirl;  increasing  a  yields  larger 
x  and  smaller  overshoot.  From  results  of  numerical  integration,  for 
finite  a,  x-  approach  finite  values  as  x^  •>  °°  for  fixed  0.  For  the  plausible 
values  ct  =  6  =  0.2,  for  x^  =  2,  5,  10,  the  corresponding  values  of  x  =  0.77, 
0.72,  0.70.  Hence,  swirl  overshoots  in  the  range  of  about  20%  seems  plaus¬ 
ible,  but  not  much  more;  estimates  that  swirl  in  the  turnaround  exceeds  the 
swirl  at  x  =  1,  z  =  h  [in  terms  of  the  definitions  of  (3.6)]  by  about  100% 
(Lewellen  1977)  are  excessive  according  to  this  analysis.  Further  results 
are  given  in  Table  1  and  graphical  presentation  is  given  in  Figures  4  and 
5. 

As  suggested  earlier,  the  periodic  nature  of  the  mathematical  solution 
is  not  observed  physically;  rather  "vortex  breakdown"  (Hall  1972;  Leibovich 
1978)  intervenes,  and  no  attempt  is  undertaken  here  to  describe,  via  the 
inviscid  mode 1  developed  in  this  section,  details  of  the  breakdown  process. 
Rather,  the  following  inferences  are  drawn  concerning  two-cell  vortical 
structure.  An  inviscid  theory  suffices  to  describe  the  nearly  horizontal 
separation  of  the  bulk  of  the  low-level  inflow  layer,  which  turns  upward 
(with  a  maximum  local  swirl  speed  probably  not  much  more  than  20%  greater 
than  the  peak  outer-vortex  value),  undergoes  a  breakdown,  and  emerges  as 
a  vertical  free  layer.  The  further  evolution  with  altitude  of  the  layer 
is  treated  in  the  next  section  by  a  model  which  retains  buoyancy,  density 
variation,  and  change  of  angular  momentum  with  height;  that  model  is  intended 
to  hold  through  a  considerable  depth  of  the  troposphere,  if  not  all  the  way 
to  the  tropopause. 

The  relatively  thin,  diffusive,  very-low-level  subregion  of  the  inflow 
layer  may  possibly  continue  inward  to  the  axis  and  fill  the  very  bottom 
of  the  central  region,  if  the  compress ionally  heated  downflux  in  the  eye 
does  not  extend  from  tropopause  all  the  way  to  surface  level.  This 
"filling  flow"  emanating  from  thin  diffusive  inflow  sublayer  may  be 
plausibly  pictured  as  a  flux  that  at  greater  height  serves  as  a  transi¬ 
tion  between  (1)  the  modest  recirculatory  motion  in  the  relatively  dry 


and  nonrotating  eye,  and  (2)  the  relatively  rapidly  rotating  and  ascending 
air  of  the  eye  wall,  which  is  characterized  by  a  near-moist-adiabatic 
locus  of  thermodynamic  states.  Indeed,  the  mass  and  momentum  interchange 
owing  to  the  shearing  interaction  along  the  fairly  well-defined  free 
boundary  between  the  eye  and  the  eye  wall  constitutes  a  possibly  fruitful 
line  of  future  inquiry  concerning  core  structure  in  a  mature  severe  vortex. 
The  entrainment  process  between  a  fast-moving  stream  and  an  adjacent, 
nearly  stagnant  stream  is  not  fully  understood,  and  existence  of  a  signi¬ 
ficant  swirling  component  of  velocity  in  the  high-speed  stream  only  further 
complicates  the  subject.  This  matter  is  taken  up  again  later  when  sugges¬ 
tions  about  future  directions  of  investigation  are  set  forth,  but  the  above 
preliminary  comments  at  this  point  seem  apropos. 
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4.  A  ONE-LAYER  TREATMENT  OF  THE  EYE  WALL 

While  some  repetition  is  inevitably  entailed,  a  self-contained 
derivation  of  a  one-layer  treatment  of  the  eye  wall  that  holds  through 
the  depth  of  the  troposphere  is  now  given.  It  is  useful  as  a  start  to 
review  some  tenets  of  the  overall  model  of  the  mature  vortical  storm. 

The  general  character  of  the  configuration  is  given  in  Figure  6. 

In  region  I  any  radial  and/or  vertical  motion  is  supposed  to  be  so  small 
that  it  is  ignorable  in  the  dynamic  balance.  The  angular  momentum  F*  is 
prescribed,  most  simply  as  a  constant,  but  probably  more  pertinently  as  a 
function  of  altitude  z*.  The  ambient  thermodynamic  state  holding  at  each 
z*  at  r*  is  taken  as  known.  Region  II  is  the  boundary  layer  in  which 
swirling  fluid  moves  radially  inward;  r*  is  the  radius  near  ground  level 
at  which  the  pressure  in  the  potential  vortex  of  region  I  is  the  same  as 
the  pressure  at  r*  =  0,  z*  =  0  when  a  fully  developed  eye  extends  from  the 
tropopause  altitude  z*  to  z*  =  0.  Region  III  contains  swirling,  updraft  fluid 
and  its  state  trajectory  is  dynamically  moist  adiabatic;  i.e.,  the  total 
enthalpy  (based  on  sea-level  ambient  conditions)  is  constant.  Region  IV 
is  nearly  stagnant  and  contains  recompressed  air  from  z*,  the  altitude  at 
which,  in  the  ambient,  the  total  enthalpy  is  the  same  as  that  at  surface 
level.  Typically,  that  ambient  enthalpy  has  a  profile  with  altitude  of 
the  character  depicted  in  Figure  7.  Recompressed  air  from  z*  may  not 
entirely  fill  region  IV  for  cases  in  which  the  insertion  of  an  eye  is 
incomplete;  the  lower  portion  of  region  IV  is  then  characterized  by 
moist-adiabatic  air,  and  the  pressure  at  r*  =  0,  z*  =  0  is  higher  than 
that  holding  for  a  completely  developed  eye.  The  locus  of  thermodynamic 
states  characterizing  the  portion  of  region  IV  filled  with  recompressed 
air  is  dry-adiabatic,  with  the  (moisture-free)  state  at  z*  serving  as  a 
reference;  the  thermodynami c  state  holding  for  the  lower  portion  (if  any) 
of  region  IV  not  filled  with  recompressed  tropopause-level  air  is  obtained 
from  hydrostatics  and  the  pressure-density  pairings  for  the  previously 
discussed  moist  adiabat,  the  pressure  and  altitude  (but  not  the  density) 
being  continuous  at  the  interface  between  dry-adiabatic  air  and  moist- 
adiabatic  air.  It  has  already  been  noted  in  previous  sections  that  the 
site  of  region  III  must  be  found  in  the  course  of  solution. 
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Finding  the  locus  of  thermodynamic  states  that  characterizes  the 

moist  adiabat  based  on  sea-level  ambient  conditions,  assigning  an  altitude 

to  the  tropopause  z*,  and  determining  the  pressure  distribution  with  height 

in  region  IV,  denoted  p*  (z*),  are  standard  steps  (Fendell  1974).  However, 

eye 

some  remarks  about  the  moist  adiabat  are  presented  in  Appendix  A. 


In  a  generalization  of  previous  notation,  the  pressure  in  region  I 
is  denoted  p^r^z*),  where  pt(r*,z*)  =  Pamb(z*)>  given,  where  "amb" 
abbreviates  "ambient"  and  the  two  subscripts  are  used  interchangeably. 
In  conventional  notation,  if  the  subscript  unity  is  discarded  for  the 
moment  for  brevity,  in  region  I, 


hence. 


9*  P*  ■--=*- 


'*2  * 
riT  P* 

r  Z* 


If  one  seeks  solution  in  the  form 


i.e.,  without  loss  of  generality, 

s*(r*)  = 


m*(z*)  = 


-f 


,*2/r*  t 

(4.1) 

l*  ; 

(4.2) 

r*2(z*)  p* 

(4.3) 

yw  ★  3  ^ 

;)  -  s*( r*)j  , 

(4.4) 

nary  derivative) 

r*2(z*)m*’(z*) ]  =  0  ; 

(4.5) 

1 

"  f 

(4.6) 

l*r*2 

dz  ’ 

(4.7) 

r*2 ( z 1 ) 


For  the  special  case  of  a  vortex  in  region  I  that  is  invariant  with 
altitude,  r*  =  r*,  const.. 
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(4.8) 


m*(z*)  =  z*/r*2; 


the  choice  of  p*(n*)  which  matches  the  ambient  atmosphere  at  r*-2  <<  1  is 
(restoring  the  subscript  unity) 


n*  =  n  * 

Ki  ^ambient 


z*  +  r*2 


/(2g*r*2)] 


(4.9) 


[In  this  approximation  the  ambient  formally  is  taken  to  hold  at  an  infinite 
radial  distance,  but  the  discrepancy  from  a  match  to  an  ambient  at  a  value 
of  (r*/r*)  of  only  ten  results  in  an  error  of  about  merely  one  percent.] 

For  the  somewhat  more  plausible  case  of  a  vortex,  the  angular  momentum 
of  which  decreases  linearly  from  a  value  of  r*  at  sea-level  to  zero  at  the 
tropopause  z*  =  z*. 


such  that  at  r*  ■*  », 


*  =  + 
2*  -  Z* 
t 


2g*r*2z£ 


(4.10) 


z*(n*  -  l) 


(4.11) 


Thus  a  solution  of  the  first-order  partial  differential  equation  p*(n*) , 
any  function  of  n*,  also  satisfies  the  boundary  condition, 

pr(r|*>  *  ^blent'2*’  for  r*  *  *• 


if  one  takes  as  the  function 


p: ^  ^  ^ambient 


rj  (Mz*  -  z*) 
2g*r*2Mz* 

:*2(mz;  -  z*) 
0  t _ _ 

2g*r*2Mz*2 


(4.12) 


Here  M  is  an  assignable  parameter  in  the  range  ^  <  M  <  1 .  Larger  values  of 
M  reduce  the  effective  dependence  of  the  angular  momentum  on  z*  and,  indeed, 
as  M  *  «  this  equation  degenerates  to  equation  (4.9). 
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*4 


More  generally, 


if  r*(z*)  =  r*[(z*  -  z*)/z*]n,  then 


(z*  -  z*) 


/  ,  r*2 

H  -  2n) _ *0 


(z*  -  z*)1_2n  2g*r*2z*2n 


i-2n 


,  n  f  j 


( z*  -  z*)  exp  ^-r*2/(2g*r*2z*) 


9  n  -  2 


(4.13) 


Finally,  if  r*(z*)  =  F*  exp(-a*z*),  then 


p*(n*)  =  p*  . 

v  J  *amb 


z*  +  (2a*)"1 in 


1  + 


rfa* 


g*r*2  exp(2a*z*) 


(4.14) 


Thus,  the  deviation  of  the  isobars  from  their  asymptotic  form  as 

horizontal  planes,  owing  to  rotation,  is  known  throughout  region  I  for 

several  stratifications  of  the  angular  momentum.  This  deviation  is  of 

significance  in  ascertaining  the  pressure  differential  at  any  altitude 

between  the  edge  of  the  eye  wall  contiguous  to  the  vortex  of  region  I 

and  the  edge  of  the  eye  wall  contiguous  to  the  eye  of  region  IV,  where 

the  eye  wall  is  envisioned  as  a  relatively  narrow  layer.  More  explicitly, 

in  region  III  the  pressure  on  the  streamline  contiguous  to  region  IV  must 

match  the  specified  p*  (z*)  in  the  eye,  and  the  pressure  on  the  stream- 

eye 

line  contiguous  to  region  I  must  match  the  pressure  in  the  potential 


vortex  p*  .  . 

Hambi ent 


j^m*-1  (m*(z*)  -  s*(r*)) 


These  statements  constitute 


the  boundary  conditions  for  the  eye-wall  calculation,  to  which  attention 
is  now  turned. 


The  symbol  R*(z*)  henceforth  denotes  the  radial  displacement  of  the 
outside  surface  of  the  "eye  wall”  with  structure;  the  symbol  h*(z*) 
denotes  the  transverse  thickness  of  the  eye  wall,  though  the  eye  wall  is 
taken  as  thin,  so  h*(z*)  <<  R*(z*).  The  velocity  component  in  a  plane 
containing  the  axis  of  symmetry  continues  to  be  denoted  by  q*,  while  the 
velocity  component  in  a  plane  perpendicular  to  the  axis  of  symmetry 
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continues  to  be  denoted  v*.  The  angular  momentum,  in  general  r*v*(r*,z*) , 
is  denoted  r*(z*);  at  r*  =  r*,  the  outer  edge  of  the  vortex,  the 
static  pressure  p*(r*,z*)  p*(r*,z*)  =  Pgmb(2*)*  a  measured  quantity, 

where  r*  »  R*(z*).  The  pressure  at  the  interface  between  the  eye  wall 
and  the  inviscid  vortex  of  region  I  is  given  by  p*[R*(z*) ,z*]  = 


p*mb{m*-1[m*(z*)  -  s*(r*(z' 


regio 

')>]}. 


where  the  formulae  for  m*(z*)  and 
The  pressure  at  the  interface  between  the  eye 

The  symbol  z*  denotes  the  altitude 
down  to  which  hydrostatic,  adiabatic  compression  (of  air  whose  with 
reference  thermodynamic  state  is  that  of  the  ambient  at  z*)  defines 


F*(z*)  are  given  above. 

and  the  eye  wall  is  given  by  p*  (z*) 

eye 


p*  (z*),  where,  again,  z*  is  the  altitude  of  the  tropopause  (the  altitude  at 
eye  t 

which  the  pressure  and  temperature  of  the  ambient  match  that  of  ground- 
level  air  which  has  passed  through  the  locus  of  thermodynamic  states 
referred  to  as  the  moist  adiabat).  For  z*  >  z*  >  Q,  P*ye(z*)  is  defined 
by  hydrostatics,  and  by  moist-adiabatic  states  such  that  the  pressure  is 
continuous  at  z*  =  z*.  If  z*  =  0,  then,  in  R*(z*)  >  r*  >  0,  the  dry- 
adiabatic  relation  holds  through  the  entire  depth  of  the  troposphere, 
defined  here  as  z*  £  z*  0. 


It  is  perhaps  also  worth  reiterating  that  standard  thermohydrostatic 
calculation  gives  P|ye(z*)  and  z\  from  specification  of  PambU*)'>  in  the 
course  of  such  a  calculation  the  triplet  ( pm0ist,Tmoist,pmoist^  has  been 
associated,  where  P*0^st  ranges  between  P*mb(0)  and  P^mb^t^’  and  the 
subscript  moist  designates  moist-adiabatic  [defined  as  the  dry  adiabatic, 
with  constant  water-vapor/dry-air  mixing  ratio,  until  saturation  and  as  the 
conventional  pseudo-adiabat  above  saturation  (see  Appendix  A)]. 


It  is  convenient  to  introduce  the  concept  of  a  series  of  calculations, 

i  =  1,2,3...  for  known  Pamb(2*)>  Tamb^z*^  pamb^Z*^  and  given  These 

calculations  are  based  on  different  completeness  of  "flushing"  of  moist- 
adiabatic  air  from  the  eye  as  denoted  by  z*^\  i  =  1,2,3...  .  For  example, 
one  could  study 


<(  i )  _  i  -  1 


*t’ 


i  =  1,2,3, 


N  +  1  , 


(4.15) 


where  i  =  1  would  constitute  an  eye  fully  flushed  of  moist  air,  i.e.,  a 
completely  developed  dry-adiabatic  eye,  and  i  =  N  +  1  would  constitute  a 
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degenerate  case  of  no  flushing,  i.e.,  no  insertion  of  dry-adiabatic  air 
in  the  core.  It  may  be  anticipated  that  the  theory  advanced  here  is 
limited  to  cases  in  which  z *  is  a  large  fraction  of  z£;  i.e.,  attention 
is  at  the  outset  confined  to  fully,  or  near  fully  developed,  eyes. 

It  is  perhaps  also  worth  a  brief  degression  to  point  out  explicitly 
that  the  use  of  hydrostatics  constitutes  a  limitation  on  what  the  present 
approach  can  ascertain  about  eye  insertion,  which  in  its  early  stages  may 
well  require  retention  of  dynamical  considerations.  A  fruitful  direction 
may  be  modification  of  conventional  plume  theory,  which  yields  the  ver¬ 
tical  profiles  of  transversely  averaged  quantities  for  buoyant  axisymmetric 
columns  over  maintained  low-level  mass,  momentum,  and  heat  sources.  The 
modification  would  entail  generalization  to  spin-up  in  swirling  environ¬ 
ments,  since  the  conventional  theory  is  for  steady  plumes  in  a  nonrotating 
environment;  what  efforts  exist  for  swirling  environments  do  not  undertake 
temporal  evolution  of  spin-up.  The  modification  would  also  entail  inclu¬ 
sion  of  provision  for  a  downdraft  core  enveloped  by  an  updraft  annulus, 
since  conventional  plume  theory  is  for  so-called  one-celled  structure  only 
(updraft  everywhere). 

It  is  useful  to  introduce  the  peak  swirl  speed  of  a  potential  vortex 
(with  a  nonrotating  core)  supportable  by  the  surface-level  pressure  deficit 
between  eye  and  ambient,  under  the  cyclostrophic  approximation: 


PaV01  * 


»e>>  1  7 


amb 


(0)  v*(l) 


(4.16) 


where  taking  the  density  fixed  at  its  ambient  value  is  adequate  for  present 
purposes,  and  the  superscript  i  is  recalled  to  denote  the  completeness  of 
eye  flushing.  For  i  =  1,  the  swirl  speed  is  maximum,  and  this  value  is 
used  for  nondimensional i zation  below. 

The  reference  value  for  the  angular  momentum  F*  is  taken  to  be 

r*  •  r*  V*(l),  (4.17) 


where  r*  is  the  (specified)  radius  at  which  the  surface  inflow  layer 
separates.  If  F*  is  held  fixed  for  all  i,  then,  if 
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(4.18) 


r*(l)\j*(l)  -  v*  ^ 

0  v  -  lo  ^  ro  Vv 


it  follows  that 


r*(1)  =  r 


rO) 


l  ’ 


>  ro 


(i) 


since  V' 


,(D 


>  v* 


(i) 


i  >  1 


(4.19) 


Also,  reference  values  for  the  streamwise  and  swirl  components  of 
velocity  in  the  eye  wall  are  taken  as  follows: 

q:(i)  =  (0.7)V*(i),  v* 

u  0 

It  may  be  noted  that 

yt + y(t> 

while  such  an  equality  holds  in  the  inviscid  portion  of  the  surface  inflow 
layer  and  hence  in  the  low-level  turn-around,  no  such  equality  is  appro¬ 
priate  for  the  eye  wall  after  the  losses  associated  with  the  breakdown. 

The  factors  (0.7)  and  (0.5)  are  rough  intuitive  guesses,  and  better 
insight  is  not  anticipated  to  be  available  in  the  near  future. 

By  conservation  of  angular  momentum  for  the  eye  wall, 

|r*^(z*)J  £v*(i)(z*)j  =  r*(i)v*(i\  (4.22) 

As  noted  earlier,  nondimensional ization  is  taken  as  follows: 


=  (0.5)V*^ 


(4.20) 


,0) 


(4.21) 


,*<*> 


r(1)(z)  ■ v(l)(z)  =1^nr 

'  0 


1 ,  =  -  2* 
V*vtj  ’  -*(  1  ) 

0 


(4.23) 


R(i)(z) 


v(0(z) 


=  I | ^ ,  where  1 ^ 


0  0  .  n  p 

:.vrv:m  -  °-5- 


(4.24) 
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By  conservation  of  mass  for  the  eye  wall, 


[°;oist(2*>]  [<!*(,><2*>]  [h*lf)<2*)_ 


=  p*^(0)a*^r*^  ^h*^1  ^ 
peye  'u'qo  ro  no 

(4.25) 


where,  for  a  thin  layer,  it  is  reasonable  to  adopt 


h*^1 ^  =  <r*,  given, 

o  o  3 


where  r*  =  300-500  miles  for  a  hurricane,  r*  =  6-8  miles  for  a  tornado, 
and  k  =  0.0025  for  both.  If  one  defines 


„  -  moist 

pmoist^  '  (l ) 


(z*) 


,(i ) 


then 


Pv  ;(0) 

eye' 


,7)  =  u(i  )/?\  _  h*(i)(z*) 

(z)  -  ^nr  *  (z)  ‘  r*(o  • 

0 


(4.26) 


where 


pmoist^ 


q(0(z) 


R(l)(z) 


h(i)(z) 


P*(.i}(0)  u*( i )  r*(i) 


3  y(i) 

L2 


<r; 


.(f)  -  [0  ,,  ""  V*' 

2  "  ^  ^TTT 

eye  o  o 


(4.27) 


p(l)(0)  r* 
(0.7)<^f-  -fij. 

-r7 


(4.28) 


Since  there  is  taken  to  be  a  constant  head  across  all  streamlines  of 
the  eye  wall,  Bernoulli’s  equation,  generalized  from  the  form  employed 
earlier  for  the  turnaround  by  the  inclusion  of  density  variation  (where 
the  moist-adiabatic  locus  relating  pressure  and  density  is  recalled  to  be 
appropriate  for  the  eye  wall)  and  of  gravitational  effects,  is  (g*  is  the 
magnitude  of  the  gravitational  accleration) 


M 

J 
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dp* 


°moist^pi ' 


(4.29) 


If  one  defines 


where 


P* 


(i) 


(z*) 


and 


P 


(i) 


(z) 


(0) 


(4.30) 


(4.31) 


P 


★ 

ref 


(4.32) 


then 


[qCl)(z)]2  *  [v(,)(z)J2  +  e3  z  ♦  j 


P«>(z> 


.-dPl  .  ■  t<i> 

pmoist(p!) 


where 


2g*r 


*».*  ( l ) 


ref 


p*^(0)/p*^(0) 

=  Keye  v  ''eye 


K’T 


/z 


(4.33) 


=  [(0.7)2  +  (0>5)2 


k(*Y 


LV 


v*' 

THT 


(4.34) 


While  the  dimensionless  form  of  the  Bernoulli  equation  just  derived 
is  the  form  useful  in  the  computations  described  below,  it  seems  worth¬ 
while  to  note  implication  derived  by  differentiation  of  the  dimensional 
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(4.35) 


(4.36) 


(4.37) 


In  general,  for  fixed  z*,  p*mb  >  P*oist  >  P*ye;  whereas  p*mb  appreciably 

exceeds  P^01-St  probably  only  for  tornado-prone  ambients  (e.g.,  cold  air 

overriding  moist  air  in  the  American  Midwest),  P*oist  appreciably  exceeds 

o*ye  except  for  a  fairly  deep  base  of  moist-adiabatic  air  in  the  very 

near-axis  region  of  the  vortex  (e.g.,  hurricanes,  which  in  fact  usually 

do  not  have  fully  dry,  fully  developed  eyes).  The  point  is  that  the 

factor  (p*^  +  p*  ,  -  2p*  .  may  be  negative,  such  that  the  two-cell 

structure  being  outlined  may  prove  compatible  with  only  certain  ambient 

stratifications;  for  if  the  just-cited  factor  is  negative,  [T*^] 

Use1/  may  turn  negative,  i.e.,  the  postulated  two-cell  structure  cannot 

T  (i)"l2 

be  realized.  Thus,  the  relation  for  t*'  '  z*  not  only  furnishes  a 
means  o'  checking  numerical  accuracy,  but  also  provides  a  hint  that  the 
present  development  may  be  able  to  distinguish  between  circumstances  com¬ 
patible  with  two-cell  structure  and  circumstances  incompatible  with  such 
structure. 
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Finally,  the  balance  of  the  pressure  gradient  across  the  eye-wall 
layer  with  accelerations  in  the  layer  yields  an  equation  that  is  a  slight 
generalization  of  one  given  earlier  for  the  turnaround: 


z  Plmb(z;R(l>)  - 


(z) 

1 2 

(Rin! 

f 

1/2 

[,<'> 

(z)' 

2  n(  i  ) 

Rzz 

2' 

3/2 

(4.40) 


Conservation  of  angular  momentum,  conservation  of  mass,  Bernoulli's 
equation,  and  the  balance  of  forces  just  stated  constitute  four  equations 
for  the  four  dependent  variables  R^(z),  v^(z),  q^(z),  and  h^(z). 

Four  dimensionless  parameters  must  be  specified:  E^,  ^3»  E4»  and  ; 
further,  the  stratification  of  a  storm-prone  (convectively  unstable) 
ambient,  e.g.,  relative  humidity  and  temperature  as  a  function  of  pressure 
through  the  depth  of  the  troposphere,  must  be  specified,  so  that  thermo¬ 
hydrostatic  calculations  giving  zt,  P  e(z),  pmoiSt^z^  as  a  function  of 
Pmoist(2)  [later  identified  with  H(Peye  +  Pamb)L  etc.,  may  be  carried 
out  Aside  from  universal  constants  such  as  g*,  specification  of  the 
four  above-cited  dimensionless  parameters  actually  entails  specification 
of  the  reference  thickness  of  the  inflow  layer,  xr*;  the  radial  distance 
from  the  axis  at  which  the  surface  inflow  layer  separates,  r*;  the  height  to 
which  moist-adiabatic  air  remains  in  the  core,  z£^ ;  the  variation  of 
the  angular  momentum  of  the  inviscid  vortex  r*  with  height,  z*,  where 

r*  =  F*  at  the  surface  z*  =  0,  and,  to  reiterate,  F*  =  r*V*^ ,  V*^ 

0  0  1 

being  the  near-surface  peak  swirl  speed  of  a  potential  vortex  with  a  non¬ 
rotating  core  (r*  >  r*  >  0)  and  a  fully  developed  eye,  for  a  given  ambient 
thermodynamic  state. 

The  set  (4.24),  ( 4 . 27 ) - ( 4 . 28 ) ,  (4.33}-(4.34) ,  and  (.4.40)  can  be 
reduced,  by  substitution  of  the  first  three  into  the  last,  to  one  second- 
order  nonlinear  ordinary  differential  equation  for  R(z);  admittedly,  the 
factor  q(z)  is  expressed  in  terms  of  an  intricate  functional  of  R(z),  but 
the  reduction  is  still  convenient  for  numerical  integration  on  a  high¬ 
speed  digital  computer.  Such  integration  requires  identification  of  two 
boundary  conditions;  these  are  taken  to  be  (for  a  vertically  rising  eye 
wall  after  breakdown,  taken  to  occur  at  so  low  an  altitude  that  z  =  0) 
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(4.41) 


Rz ( 0 )  =  0,  R(0)  chosen  such  that  Rzz(°)  =  0  in  (4.40). 


Conceptually,  this  initial-value  problem  may  then  be  marched  forward  in  z 
from  z  =  0  to  z  =  2t  (or  even  further  if  extrapolation  of  tabulated  data 
is  accepted).  In  practice,  singular  behavior  virtually  always  intervenes, 
such  that  some  physical  condition  of  the  formulation  is  violated,  and  fur¬ 
ther  calculation  is  not  justified.  For  example,  the  layer  is  taken  thin, 
but  if,  at  any  height  z,  the  thickness  h  becomes  comparable  to  the  radial 
position  R,  the  approximations  are  inappropriate. 

It  may  be  remarked,  with  the  completion  of  the  formulation  of  the  one- 
layer  theory  for  the  eye  wall,  that  proceeding  to  more  highly  resolved 
transverse  structure  is  not  a  trivial  extention.  For  example,  proceeding 
to  a  two-layer  model  entails  distinguishing  the  significant  transverse 
variation  of  the  pressure  field  in  the  Bernoulli  equation  pertinent  to 
each  sublayer.  But  the  total  head,  the  swirl,  the  streamwise  velocity 
component,  and  the  gravitational  term  could  be  the  same  for  each  of  two 
comparably  thick  sublayers  constituting  the  relatively  thin  eye  wall,  so 
such  distinguishing  requires  subtlety.  Multiple-subdivision  analysis  is  not 
pursued  here  and  attention  is  confined  to  the  one-layer  formulation. 


If,  in  the  one-layer  formulation,  one  has  the  same  thermodynamic 
stratification  and  angular-momentum  stratification  of  the  ambient,  the 
same  degree  of  flushing  of  the  eye,  and  the  same  ratio  of  inflow-layer 
thickness  to  inflow-layer-separation  distance  for  a  fully  developed  eye 
(h*^/r*^),  then  the  only  distinction  in  the  dimensionless  formulation 
between  two  cases  (one  of  which  conceivably  could  be  a  hurricane  case  and 


the  other  a  tornado  case,  within  the  above  constraints)  is  the  dimension¬ 
less  parameter  Z3  £=  (2g*r*^ 1 ^)/[ V*^ 1 ^  J,  in  that  the  radial  scale  to 
peak  low-level  wind  position  for  a  fully  developed  eye,  r*^  =  r*, 
might  be  quite  different  in  the  two  cases.  The  parameter  l  enters  in 


the  gravitational  term  of  the  Bernoulli  equation. 


30 


5.  PRELIMINARY  RESULTS 


Results  of  the  numerical  computations  are  presented  in  Figures  8-10. 
Typical  thermodynamic  ambient  data  is  analyzed  for  a  hurricane  environment 
[specifically,  a  typical  West  Indies  environment  for  September  (Jordan 
1948)]  and  a  tornado  environment  (unpublished  sounding  taken  in  Jackson,  MS, 
April  1978,  in  spatial  and  temporal  proximity  of  a  moderate-intensity 
tornado).  For  each  ambient  profile,  the  thermodynamic  variables  (as 
functions  of  z*)  characterizing  the  moist  adiabat  and  eye  region  are 
calculated  as  described  in  Appendix  A  (see  also  Fendell  1974),  as  a  pre¬ 
requisite  to  the  eye-wall  calculation  described  in  the  previous  section. 

For  both  hurricane  and  tornado,  under  various  adopted  stratifications 
for  the  ambient  circulation,  the  initial  value  of  R  (the  outer  edge  of 
the  eye  wall  at  ground  level,  z  =  0)  increases  monotonical ly  with  an 
increase  in  the  altitude  to  which  the  axial  region  (region  IV,  Figure  6) 
is  filled  with  moist  air.  This  result  may  be  anticipated  as  follows: 
an  increase  of  the  height  to  which  moist  air  persists  in  the  nominally 
dry-adiabat  eye  region  reduces  the  surface-level  pressure  deficit  between 
the  axis  and  the  ambient.  Reduction  in  pressure  deficit  results  in  a 
decrease  of  the  peak  swirl  speed  of  the  main  vortex  under  the  cyclostrophic 
approximation  (4.16).  It  then  follows,  from  the  conservation  of  angular 
momentum,  (4.24),  that  R  increases. 

For  the  hurricane  ambient,  the  observed  shape  of  the  eye  wall  is 
s  jnificantly  altered  by  the  stratification  (i.e.,  the  change  with 
altitude)  of  the  angular  momentum  r  at  the  periphery  of  the  vortex 
(Figure  8).  When  this  angular  momentum  is  taken  invariant  with  altitude 
[M  -  in  (4.12)3,  the  eye  wall  shape  is  seen  to  be  nearly  vertical  at 
low  altitude,  followed  at  high  altitude  by  a  pronounced  curvature  radially 
outward  into  the  ambient  region  I;  the  eye-wall  position  R  increases  with 
altitude  z  for  higher  values  of  z.  As  the  depth  of  the  moist-air  layer 
at  the  base  of  the  eye  region  increases  [i.e.,  increasing  i  in  (4.15)]  the 
vertical  portion  of  the  eye-wall  locus  extends  further  up  in  the  tropo¬ 
sphere.  For  zj  =  0.5  z*  [see  (4.15)],  the  locus  terminates  horizontally 
(R  increasing  at  fixed  z)  at  an  altitude  about  two-thirds  of  the  distance 
to  the  tropopause.  Here  zj  is  that  cutoff  height  down  to  which  the  dry 
adiabat  is  extended  from  the  tropopause,  such  that  moist-adiabatic  air 
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lies  in  zj  >  z*  >  0;  z£  is  the  top  of  the  tropopause,  which,  for  the 
hurricane  data,  is  computed  to  be  approximately  45,000  feet.  As  more 
adiabatic  air  is  flushed  from  the  eye  (i.e.,  for  zj  decreasing),  the 
termination  altitude  also  decreases  until,  for  an  eye  completely 
flushed  of  moist  air  so  that  zj  =  0,  the  termination  altitude  is  only 
about  one-fifth  z*. 

If  the  angular  momentum  in  the  vortex  is  functionally  prescribed 
to  decrease  to  zero  at  z*,  a  marked  change  in  the  computed  eye-wall 
shape  is  noted.  In  general,  the  faster  this  decrease  of  r  with 
altitude  (i.e.,  the  smaller  M  in  (4.12)],  the  stronger  is  the  tendency 
for  the  eye  wall  to  penetrate  radially  toward  the  eye  center.  For 
M  =  1,  the  "eye  wall  profiles"  actually  turn  horizontally  toward  the 
vortex  axis. 

For  the  more  moderate  value  M  =  2,  elongated  ”S”  shaped  profiles 
are  produced.  These  eye-wall  loci  penetrate  higher  into  the  tropopause 
than  both  those  for  smaller  values  (M  =  1)  and  larger  values  (M  -*•“). 
Indeed,  for  z*  =  0.5  z*  [i.e.,  for  i  =  6  in  (4.15)],  the  eye  wall  reaches 
the  top  of  the  tropopause  (z  =  45,000  ft).  As  discussed  in  the  next 
section,  future  work  will  involve  seeking  to  identify  those  angular- 
momentum  stratifications  which  can  yield  a  vertical  eye-wall  locus 
through  the  entire  extent  of  the  troposphere. 

In  Figure  9,  the  altitudinal  variation  of  the  eye-wall  thickness  h, 
the  velocity  component  in  a  plane  containing  the  axis  of  symmetry  q,  and 
the  velocity  component  (swirl)  in  a  plane  perpendicular  to  the  axis  of 
symmetry  v,  are  presented.  The  case  examined  is  "unstratified"  (i.e., 

M  -*•  °°) ,  with  the  eye  consisting  of  moist  air  from  the  ground  up  to 
zj  =  0.5  z*  [i.e.,  i  =  6  in  (4.15)].  The  eye-wall  thickness  is  found 
to  increase  monotonical ly  with  height.  The  theory  is  based  on  a  thin 
eye  wall,  i.e.,  h  «  R  at  all  z,  and  indeed,  h  <_  0.22  throughout,  while 
R  >  1.7  everywhere.  The  swirl  velocity  v  is  nearly  constant  at  low 
altitude,  then  decreases  at  higher  altitudes.  This  behavior  follows 
directly  from  the  radially  outward  bending  of  the  eye-wall  location  R, 
since  this  particular  case  entails  angular  momentum  r  =  R  v  =  const. 
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Results  for  the  tornado  ambient,  Figure  10,  exhibit  a  characteris¬ 
tically  different  solution.  Here,  for  each  of  the  angular-momentum- 
stratification  parameters  examined  (M  =  1,  2,  °°) ,  the  completely  dry 
eye  (i  =1)  represents  a  "ceiling"  solution  for  the  eye-wall  shape. 

That  is,  the  eye-wall  solution  locus  monotonically  increases  radially 
outward  with  increasing  altitude  and  represents  an  upper  bound  such  that, 
for  a  given  M,  resultant  eye-wall  loci  for  increasingly  moist  eyes  do  not 
penetrate  much  above  this  "ceiling"  locus.  As  a  deeper  layer  of  moist 
air  is  taken  to  exist  at  the  base  of  the  eye  region,  i.e.,  as  i  increases 
from  unity  to  larger  values,  the  resultant  eye  wall  loci  rise  from  their 
surface-level  values  until  they  meet  the  "celling"  solution  at  approxi¬ 
mately  x*.  They  then  roughly  follow  this  ceiling  solution  radially 
outward,  such  that  R  increases  to  larger  values. 

Like  the  hurricane,  the  radial  displacement  R*  of  the  outside 
surface  of  the  eye  wall  at  the  ground  z*  =  0  increases  monotonically  as 

additional  moist-adiabatic  air  is  added  to  the  base  of  the  dry-adiabat 

locus  near  the  axis.  Also,  as  the  angular  momentum  of  the  peripheral 

vortex  decreases  more  rapidly  with  height,  i.e.,  as  M  -*■  1  ,  the  eye  wall 

shape  "bends"  more  acutely  toward  the  axis  of  rotation  before  inter¬ 
secting  the  "ceiling"  solution. 
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6.  FUTURE  DIRECTIONS 


An  inviscid  axisymmetric  nonentraining  nonheat-conducting  analysis 
of  the  eye  wall  in  a  two-cell  vortex  has  been  set  forth.  The  flow  exterior 
to  the  eye  wall  is  a  potential  vortex,  and  the  flow  in  the  compressionally 
heated  eye  (interior  to  the  eye  wall)  is  static.  Preliminary  results 
suggest  a  key  role  for  the  variation  of  the  angular  momentum  of  the  outer 
vortex  r  with  altitude  z,  for  perhaps  only  with  a  limited  class  of  such 
stratifications  can  a  two-cell  structure  be  realized  for  a  given  convect- 
ively  unstable  ambient,  with  a  fully  developed  (or  fairly  well  developed) 
eye  and  reasonable  vortex  dimensions.  Very  preliminary  results  suggest 
that  a  fairly  vertical  eye  wall  through  an  appreciable  fraction  of  the 
troposphere  is  not  achieved  for  outer-vortex  angular  momentum  invariant 
with  altitude.  In  fact,  a  useful  exercise  planned  for  the  future  is  to 
take  R,  the  eye-wall  position  relative  to  the  axis  of  symmetry,  fixed  for 
all  z  at  its  low-level  value  for  altitude  z  =  0  and  to  use  the  resulting  four 
previously-derived  coupled  nonlinear  (now  entirely  algebraic)  equations  to 
obtain  the  eye-wall  thickness  h(z),  its  streamwise  speed,  q(z),  its  swirl 
speed  v(z),  and  F(z);  hopefully  for  pertinent  parametric  input,  this  inverse 
problem  [for  which  the  conventional  roles  of  F(z)  and  R(z)  are  reversed, 
in  that  now  r(z)  becomes  an  output  and  R(z)  an  input]  has  a  solution  over 

an  extended  range  of  altitude,  say  0,8  >  z  >  0.  A  more  difficult 

exercise  would  be  to  ascertain  the  T(z)  for  which  the  difference 

[R(z)  -  h(z)],  characteristic  of  the  locus  of  the  eye/eye-wall  interface, 
is  constant  in  z  over  much  of  the  range  zt  >  z  >  0.  The  statement  possible 
at  this  early  stage  of  investigation  is  that,  for  specified  input  (vortex 
parameters,  ambient  atmospheric  sounding  of  the  troposphere),  the  measurable 
quantity  F(z)  may  be  a  means  of  ascertaining  whether  a  one-cell  moderate 
vortex  can  become  a  two-cell,  intense  vortex;  in  fact,  no  other  candidate 
for  a  structure-differentiating  observable  seems  readily  identifiable. 

Finally,  some  brief  comments  are  appended  concerning  a  longer-range 
plan  of  investigation,  which  incorporates  dynamics  in  the  eye,  since  the 
work  discussed  thus  far  is  limited  to  a  static  description  of  the  eye. 

It  should  suffice  to  state  but  once  the  obvious  proviso  that  the  direction 
of  sequential  investigations  are  heavily  contingent  on  results  of  earlier 
investigations.  It  has  already  been  mentioned  that,  upon  adoption  of  a 


34 


family  of  very  simple  models  of  entrainment,  one  should  investigate 
whether  the  steady  supply  of  very-low-momentum  fluid,  entering  the  eye 
region  (with  sea-level  enthalpy)  from  the  much  thinner,  frictional  sub¬ 
layer  (under  the  more  energetic,  effectively  nondiffusive  outer  portion 
of  the  inflow  layer),  is  carried  away  by  the  entrainment  process.  If 
this  is  possible  at  the  postulated  eye  pressure,  then  a  geometry  consis¬ 
tent  with  the  postulated  eye  pressure  arises. 

Accordingly,  one  could  adopt  an  initial  condition  in  which  the  eye 
pressure  is  higher  than  that  of  the  steady  flow  which  emerges  from  the 
studies  of  the  last  paragraph,  for  one  of  the  entrainment  models.  With 
that  same  entrainment  model,  with  a  quasi-steady  treatment  of  the  eye, 
and  with  no  entrainment  on  the  outside  of  the  rising  air  in  the  (eye-wall) 
annulus,  one  could  try  to  formulate  and  solve  a  one-spatial-coordinate, 
time-dependent  model  of  the  motion  leading  to  time-description  of  eye 
development.  By  one-spatial -coordinate  modeling,  reference  is  being 
made  to  the  same  kind  of  transversely-averaged  analysis  as  that  executed 
earlier  in  this  report.  However,  in  the  present  context,  the  outer 
potential -vortex  mass  can  migrate,  inertially,  and  the  proportions  of 
"tropopause  air"  and  sea-level  air  in  the  eye  can  change.  The  inertial 
migration  refers  to  an  oscillation  established  during  the  later  stages 
of  tropical-cyclone  intensification,  as  radially-inward-moving  particles 
overshoot  their  radial  position  of  equilibrium  (in  which  centrifugal 
acceleration  is  made  compatible  with  moist-adiabatic-ascent-associated 
low-level  pressure  reduction  from  ambient  value).  The  reconciliation 
during  spin-up  of  the  radial  conservation  of  momentum  (gradient-wind-type 
equation)  and  the  axial  conservation  of  momentum  (hydrostatics )  furnishes 
constraints  discussed  elsewhere  by  Carrier  (1971a). 

It  is  worth  reiterating  that  experimental  guidance  for  the  entrain¬ 
ment  models  to  be  adopted  is  sorely  needed,  since  in  the  present  instance 
the  occurrence  of  appreciable  swirl  in  the  fast-moving  stream  alters  the 
situation  from  the  more  commonly  encountered  case,  in  which  the  fast- 
moving  stream  (contiguous  to  the  slower-moving  stream)  is  not  rotating. 

Only  after  the  core  region  of  intense  atmospheric  vortex  is  more  in 
hand  is  examination  of  the  upper-level  outflow  of  mass,  momentum,  and 
energy  from  the  typhoon,  and  its  compatibility  with  the  synoptic  situation 
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in  the  far  field,  intended.  The  pressure  deficit  associated  with  the  vor¬ 
tex  cannot  be  maintained  without  such  compatibility,  but  what  the  pressure 
deficit  in  fact  may  be  cannot  be  ascertained  without  a  better  comprehension 
of  the  one-cell/two-cel 1  transition  in  a  tropical  cyclone. 

Finally,  aside  from  matters  of  horizontal  scale,  of  lifespan,  and  of 
enthalpy  transfer  at  the  surface,  a  particularly  relevant  distinction 
between  the  exceptional  long-lived,  long-path,  high-swirl-speed  tornado 
spawned  by  a  rotating  thunderstorm  and  the  typhoon  is  the  speed  of  trans¬ 
lation  of  the  two  atmospheric  systems.  Hurricanes  drift  westward  in  the 
trades  at  about  15-20  mph ,  and,  even  in  recurvature  around  persistent  sub¬ 
tropical  anticyclones,  rarely  translate  much  in  excess  of  35  mph.  In 
contrast  tornadoes  often  translate  in  the  vicinity  of  70  mph,  and  entail 
systems  with  strong  vertical  wind  shear;  i.e.,  the  advective  component 
of  velocity  in  the  vicinity  of  the  midlatitudinal  supercell  changes 
magnitude  and  direction  with  altitude,  whereas  the  absence  of  such 
"ventilation"  is  virtually  a  prerequisite  for  the  intensification  of 
tropical  cyclones.  The  pertinence  of  the  four-part  structure  of  an 
intense  atmospheric  vortex  presented  for  tropical  cyclones  may  not  suffice 
for  the  tornado  cyclone;  especially,  the  suitability  of  regarding  the 
system  as  closed  and  containing  its  own  “captured"  airmass  (throughput 
supply)  of  convectively  unstable  air  (the  large  region  I  of  the  adopted 
four-part  structure)  is  unclear. 
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TABLE  1.  NUMERICAL  RESULTS  FOR  EXTREMA  OF  THE  TURNAROUND, 
FROM  INTEGRATION  OF  THE  INITIAL-VALUE  PROBLEM  (3.8)-(3.1G) 


a 

6 

a/S 

rioijr; 

<">,-0 

‘r>c,.c 

<‘\.o 

0.0 

0.1 

0.0 

2.0 

0.791 

0.789 

1.264 

1.266 

0.0 

0.1 

0.0 

5.0 

0.754 

0.751 

1.326 

1.329 

0.0 

0.1 

0.0 

10.0 

0.742 

0.738 

1.348 

1.352 

0.0 

0.1 

0.1 

100.0 

0.732 

0.727 

1.366 

1.372 

0.0 

0.2 

0.0 

2.0 

0.718 

0.714 

1.392 

1.397 

0.0 

0.2 

0.0 

5.0 

0.871 

0.664 

1.489 

1.497 

0.0 

0.2 

0.0 

10.0 

0.656 

0.647 

1.524 

1.533 

0.0 

0.2 

0.0 

100.0 

0.644 

0.632 

1.553 

1.565 

0.1 

0.1 

1.0 

2.0 

0.820 

0.789 

1.295 

1.266 

0.1 

0.1 

1.0 

5.0 

0.777 

0.751 

1.348 

1.329 

0.1 

0.1 

1.0 

10.0 

0.762 

0.738 

1.366 

1.352 

0.1 

0.1 

1.0 

100.0 

0.750 

0.727 

1.381 

1.372 

0.1 

0.2 

O.S 

2.Q» 

0.748 

0.714 

1.422 

1.397 

0.1 

0.2 

0.5 

5.0 

0.696 

0.664 

1.510 

1.497 

0.1 

0.2 

O.S 

10.0 

0.677 

0.647 

1.541 

1.533 

0.1 

0.2 

0.5 

100.0 

0.663 

0.632 

1.567 

1.566 

0.1 

o.s 

2.0 

2.0 

0.876 

0.846 

1.212 

1.181 

O.J 

0.5 

2.0 

s.o 

0.841 

0.818 

1.244 

1.223 

0.1 

0.5 

2.0 

10.0 

0.828 

0. 808 

1.254 

1.237 

0.1 

0.5 

2.0 

100.0 

0.818 

0.800 

1.264 

1.250 

0.2 

0.1 

2.0 

2.0 

0.845 

0.789 

1.325 

1.266 

0.2 

0.1 

2.0 

5.0 

0.796 

0.751 

1.369 

1.329 

0.2 

0.1 

2.0 

10.0 

0.779 

0.738 

1.384 

1.352 

0.2 

0.1 

2.0 
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0.765 

0.727 

1.396 

1.372 

0.2 

0.2 

1.0 

2.0 

0.774 

0.714 

1.452 

1.397 

0.2 

0.2 

1.0 

5.0 

0.716 

0.664 

1.530 

1.497 

0.2 

0.2 

1.0 

10.0 

0.695 

0.647 

1.558 

1.533 

0.2 

0.2 

1.0 

100.0 

0.678 

0.632 

1.581 

1.566 

0.2 

0.4 

0.5 

2.0 

0.685 

0.615 

1.652 

1.606 

0.2 

0.4 

0.5 

5.0 

0.617 

0,550 

1.789 

1.772 

0.2 

0.4 

0.5 

10.0 

0.594 

0.528 

1.838 

1.832 

0.2 

0.4 

0.5 

100.0 

0.578 

0.509 

1.877 

1.887 

Results  for  (x~)a=0  are  from  (3.17).  Also,  x  and  =  TOO  are 
within  2%  of  values  for  x'  -*■  for  a,8  studied. 


Figure  1 .  Schematic,  not  to  scale,  of  an  inviscid  model  of  the  turnaround 
region  of  a  severe  vertical  quasisteady  axisymmetric  vortex.  The  "eye", 
region  IV,  is  isobaric  at  pressure  p*,  over  the  vertical  extent  of  interest 
here,  where  p*  is  also  the  pressure  at  r*  =  r*,  z*  =  h*  (because  the  pres¬ 
sure  above  hydrostatic  is  approximately  invariant  across  the  surface  inflow 
layer  II).  The  "eye  wall",  region  III,  is  demarcated  by  two  streamsurfaces , 
i|i*(r*,z*)  =  const,  (the  position  of  each  to  be  determined),  encompassing 
the  mass  efflux  from  region  II.  It  is  recalled  that  a  potential  vortex 
holds  in  region  I. 
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Figure  2.  Schematic  of  the  location  of  an  inviscid  "eye  wall  without 
structure"  demarcating  the  interface  between  an  isobaric  nonswirling 
"eye"  at  pressure  p*  and  a  potential  vortex.  The  sheet  representing 
the  "eye  wall"  has  thickness  h*  -*■  0;  in  turnaround  region,  its  displace¬ 
ment  from  the  axis,  as  a  function  of  height  above  the  ground  plane,  is 
denoted  r*( z*) . 
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Figure  3.  Phase-plane  properties  of  equation  (3.12),  where  P  =  (dx/dz), 
with  x  the  dependent  variable  (x  >  0)  and  z  the  independent  variable. 
Isoclines  of  zero  and  infinite  slope  are  noted.  At  x  =  1 ,  the  distance 
from  the  axis  of  symmetry  at  which  the  surface  inflow  layer  separates, 
a  finite  positive  slope  is  adopted.  For  a  >  0,  x*  >  1;  for  a  s  0,  x*  =  1. 
The  sketched  trajectory  (solution  curve)  is  a  limit  cycle  (closed  curve 
indicative  of  periodic  behavior).  The  periodicity  is  not  of  physical 
interest. 
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Figure  4.  The  solution  to  the  boundary-value  problem  posed  by  equations  (3.9),  (3.10), 
and  (3.12),  for  =  2;  the  solution  may  be  completed  by  symmetry  considerations  to 
constitute  a  full  cycl 
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Figure  6.  Another  schematic  of  the  postulated  four-part  model  of  the 

structure  of  a  mature  intense  hurricane,  of  radial  extent  r*  and  axial 

extent  z*  (the  tropopause).  The  potential  vortex  in  region  I  is  here 

ascribed  axial  variation,  such  V  ?t  the  associated  pressure  field  is 

denoted  p*(r*,z*),  where  p*(r*,z*,  =  9iven-  The  low-level 

swirling  influx  (region  II)  grows  in  thickness,  and  separates  to  form 

the  "eye  wall"  (region  III)  of  thickness  h*;  the  flow  speed  independent 

of  swirling  is  denoted  q*.  The  pressure  variation  in  the  "eye",  region 

IV,  is  entirely  axial  and  is  denoted  p*  (z*);  for  the  idealization  of 

eye 

a  completely  dry,  cloud-free,  fully  developed  eye,  this  variation  is 
determined  by  dry-adiabatic  compression  of  tropopause-level  air. 
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Figure  7.  A  typical  profile  of  the  total  static  temperature  (H*/c*)  over 
the  ocean  in  the  equatorial  trough  and  in  the  subtropics  (1200  n  mi  from 
the  equatorial  trough),  as  a  function  of  pressure  p*  and  altitude  above 
sea-level  z*  (Palmer  and  Newton  1970,  p.  439,  575).  The  loci  of  thermo¬ 
dynamic  states  for  a  typical,  relatively  slowly  rising  cumulus  (cu)  and 
also  for  the  undiluted  core  of  a  nonentraining,  relatively  rapidly  rising 
cumulonimbus  (cb)  are  also  noted;  the  locus  (cb)  is  that  appropriate  for 
moist-adiabatic  ascent  of  sea-level  air.  The  dashed  curve  is  the  typical 
ambient  referred  to  in  the  text:  a  well-mixed  lowest  layer,  a  midtropo- 
spheric  minimum,  and  a  recovery  of  sea-level  value  at  altitude  z*  (defined 
here  to  be  the  tropopause,  and  typically  about  50,000  ft).  The  quantity 
H*  is  the  sum  of  the  static  enthalpy,  the  latent  heat  equivalent  of  the 
water  vapor  present,  and  the  gravitational  potential  energy;  the  quantity 
c*  is  the  specific  heat  at  constant  pressure,  which  may  be  taken  to  be 
effectively  that  of  dry  air. 


44 


* 


* 


Figure  8.  The  outer  edge  of  the  eye  wall,  R,  measured  from  the  center 
of  the  potential  vortex,  is  presented  as  a  function  of  altitude,  z,  for 
a  hurricane  environment.  Here  R  and  z  are  nondimensional i zed  by  r*  = 
50,000  ft.  Cases  are  presented  for  three  different  values  of  the  angular 
momentum  "stratification"  parameter  M  in  (4.12),  namely  M  =  1,  2,  and  <*>. 
For  each  case,  the  eye  wall  locus  closest  to  the  axis  (i.e.,  the  locus 
with  smallest  R  at  z  =  0),  represents  an  eye  region  completely  flushed 
of  moist  air  [i.e.,  i  =  1  in  (4.15)3.  Successive  loci  at  larger  radial 
positions  represent  eyes  increasingly  filled  with  moist  air,  specifically 
unity  increases  of  i  in  (4.15). 
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Figure  9.  For  a  hurricane  environment,  the  altitudinal  variation  of  the 

eye  wall  thickness  h,  the  velocity  component  in  a  plane  containing  the 

axis  of  symmetry  q,  and  the  velocity  component  in  a  plane  perpendicular 

to  the  axis  of  symmetry  v,  are  presented.  The  case  examined  is 

"unstratified"  (M  -*■  °°) ,  with  the  eye  consisting  of  moist  air  from  the 

ground  up  to  z*  =  0.5  z*  [i.e.,  i  =  6  in  (4.15)],  where  z*  is  the  top 

of  the  tropopause.  Here  h  is  nondimensional ized  by  r*  =  50,000  ft, 

(  1  )  1 

while  q  and  v  are  nondimensionalized  by  V*  =  422.9  ft/sec. 


I 
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Figure  10.  In  a  similar  manner  to  Figure  8,  but  for  a  tornado  environ 
ment,  the  radial  edge  of  the  eye  wall  R  is  presented  as  a  function  of 
altitude  z.  Here  R  and  z  are  nondimens ional i zed  by  r?  =  500  ft. 
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APPENDIX  A.  A  DERIVATION  OF  THE  MOIST  ADIABAT 


Consider  a  volume  V  with  surface  S  and  surface  element  dS ;  the  outward 
normal  vector  to  dS  is  denoted  n.  and  fluid  velocity  vector  is  denoted  u^ . 
Super  asterisks,  used  in  the  main  text  to  denote  dimensional  quantities, 
are  dispensed  with  for  this  appendix. 

For  conservation  of  mass,  consider  air  first.  The  flux  of  dry  air 
through  dS  is  in  Cartesian  notation,  p  u.  n^  dS,  where  p  is  air  density. 
Hence,  in  the  steady  state, 


p  ui  ni  dS  =  0  =>  (p  u.)  . 


=  0, 


( A.  1 ) 


where  comma  denotes  partial  differentiation.  Water  vapor  has  the  same 
speed,  but  density  a.  When  a  <  Z(T),  where  £(T),  the  saturation  vapor 
density,  is  a  known  function  of  the  temperature  T, 


(a  u.)  ,  =  0. 

'  »  1 

If  the  mixing  ratio  for  water  vapor  Y  h  a/p,  then 


P 


Ui 


(A. 2) 


When  T  is  low  enough. 


a  =  £(T) .  (A. 3) 

For  conservation  of  momentum,  the  momentum  flux  associated  with  u. 
is  (p  +  a)u.  n.  u.  dS,  and  that  associated  with  the  partial  pressure  of 

J  J  1 

air  p  and  with  the  partial  pressure  of  vapor  a  is  (a  +  p)n^  dS.  The  pro¬ 
duction  of  momentum  per  unit  volume  associated  with  the  gravitational 
body  force  (g  is  the  magnitude  of  the  acceleration  of  gravity)  is 

-g  z  .(p  +  a),  where  the  spatial  coordinate  z  lies  parallel  to  the  body 
*  * 

force.  The  loss  of  momentum  per  unit  volume  associated  with  the  dis¬ 
appearing  fluid  is  -(a  u.)  .  u. .  Thus,  using  (A.l), 

J  J  J  ' 
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n  .  u .  dS  + 
J  i 


(a  +  p)n..  dS  =  -g 


becomes 


u .)  .  u.  dV 

j  ,j  i 


(A. 4) 


(p  +  c)u.  u.  +  (a  +  p)  .  +  (p  +  a)g  z  .  =  0.  (A. 5) 

J  1  \  >  1  »  * 

*  J 

Incidentally,  if  a  momentum  equation  had  been  written  for  each  constituent, 
dry  air  and  water  vapor,  interparticle  drag  of  the  form  D  =  A(u^  -  v^) 
would  have  been  introduced,  where  is  air  velocity  and  v^  is  vapor  velo¬ 
city,  and  here  A  +  “  so  v.  =  u^ . 

The  equations  of  state  are 


P  =  P  R  T, 

(A. 6) 

a  »  o  R*  T, 

(A. 7) 

where  the  gas  constant  for  vapor  R'  is  related  to  the  gas  constant  for 

vapor  by  R1  =  R/(mv/ma),  where  my  is  the  molecular  weight  of  water  vapor 

and  m  is  the  molecular  weight  of  dry  air. 
d 

The  nonmechanical  energy  carried  by  the  vapor  can  be  written  as 
| L0  +  c^  t),  where  L0  is  the  (extrapolated)  value  of  the  heat  of  condensa¬ 
tion  at  T  =  0,  and  cv  is  the  specific  heat  at  constant  volume  of  the 
vapor  v.  With  this  description,  the  net  heat  released  by  condensation 
at  temperature  T  would  be 


L  =  L0  +  (cj  -  c)T ,  (A. 8) 

where  c  is  the  specific  heat  of  liquid  water.  The  flux  of  energy  across 
dS  via  u^  is 
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ni  dS, 

(A. 9) 

where  cv  is  the  specific  heat  at  constant  volume  of  dry  air.  The  stream- 
wise  velocity  component  is  denoted  q  and  the  perpendicular  velocity  com¬ 
ponent  (swirl)  is  denoted  v.  Below,  Q2  h  (v2  +  q2 ) .  The  work  done  by 
pressure  forces  across  dS  is  (a  +  pju^  n^  dS,  and  the  flux  of  potential 
energy  is  (p  +  j)g  z  n^  u..  dS .  The  rate  at  which  energy  goes  into  a  sink 
(i.e.,  into  liquid,  which,  by  conventional  approximation  in  treating  the 
moist  adiabat,  just  vanishes)  is  the  product  of  two  factors:  -(a  u..) 
the  rate  of  loss  of  mass  per  unit  volume,  and  jg  z  +  jjq2  +  v2 ) /2 J  +  c  t|. 
It  follows  that 


“f  (c 


T  + 


) t0  “i  (‘ 


T  + 


+  v" 


+  L 


I  P  u^cv  T  +  ^  +  a  u.  ^cJJ  T  +  L0  +  ^  +  (p  +  o)Uj  g  z  +  (a  +  p)Uj 
L  .  »J 

=  (a  u^)  .|c  T  +  |-  +  g  zj  ,  (A. 10) 

Hence, 

s  "J  (9  2  *  *  (a  "J>.J  [(<  -  C)T  +  Lo] 

*  °  “j(<  uj(c»  t)j 

+  (u.p.  +  pu.  )  +  (  u .  a  .  +  a  u .  1=0, 

V  j  j,j/  V  j  »j  Jj / 

(A. 11) 

where  (a  u.)  .  =  p  u .  Y  . ,  a  u .  =  p  Y  u.,  and  the  last  line  in  (A. 11)  is 


P 


(A. 12) 
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Thus,  since  (cv  +  R)  =  cp  and  (c^  +  R')  =  c^. 


u.(l  +  Y)  (g  z  +  +  p  u,(c  T)  ,  +  p  u.  Y  .(L  - 

J  \  C  /  A  JP>J  J»JO 


+  P  U.(Y  c;  T)  .  =  0, 

J  r  >  J 


(A. 13) 


where  (not  c^)  appears  because  work  is  done  on  condensing  gas  and  the 
enthalpy  release  is  ^  -  (c  -  c^lj.  But 

1  To  1  (p  +  a)  . 

P  u.  (1  +  Y)  [(Q2/2)  +  g  zj^  +  - — =  0,  (A. 14) 


so,  if  s  is  distance  along  a  streamline,  and  since 


« &  ■  “j  s:  • 

J 


(a  +  p) 


(lo-cT)-(cpt  +  Ycp  T),sa0- 


(A. 15) 


(A. 16) 


If  it  is  recalled  that 


a  =  Z(T)  R'  T,  p  =  p/R  T,  Y  =  I(T)  R  T/p, 


(A. 17) 


the  last  equation  becomes,  if  p  =  (p  +  a)  and  a  =  (mv/ma). 


_ii_ + rL .  (c .  cv)ii  s-i 

P  -  Z(T)  L  0  P  J  [p  - 

O  Z(T)  cjl  r  ~T 

5  -  mi  L  0  p  J 


o  £(T)  c;  r  v  1  0  P 

c  +  - £  +  L  -  (c  -  c^)T  - 

P  P  -  S(T)  L  0  P  J  [p  -  E(T)]2 


,  (A. 18) 


where  E'(T)  denotes  the  first  derivative,  and  all  specific  heats  have 
been  taken  as  constant.  In  some  more  approximate  derivations,  p  =  p, 
(Cp  +  Y  op  =  cp,  and  cp  =  c^  in  the  first  term  of  the  denominator. 
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